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ABSTRACT 

Despite decades-old predictions of the expected presence of dozens of cat- 
aclysmic variables in the cores of globular clusters, the number of irrefutable, 
out-bursting candidates is still barely a handful. Using multi-wavelength, multi- 
epoch HST images we have produced outburst and quiescence light curves for 
the recently discovered large amplitude variable CVl in the core of the post core- 
collapse globular cluster M15. The light curves and blue colors show that the 
object is a bona fide dwarf nova, with absolute magnitude at maximum light 
rivaling that of the most luminous known dwarf novae. 

Subject headings: Stars: dwarf novae, cataclysmic variables — globular clusters: 
individual (M15) 
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1. Introduction 

The realization that X-ray binaries are 100-1000 times overabundant in Galactic 
globular clusters relative to the rest of the Milky Way (Clark 1975; Katz 1975) was crucial 
in understanding that intimate and tightly coupled relationships exist between stellar 
evolution and stellar dynamics. Among the predictions of this theory is the expectation 
that many white dwarf/red dwarf binaries should form via two-body tidal capture (Fabian 
et al. 1975) and/or three-body exchange captures (Hut 1983; Hut & Paczynski 1984). Some 
of these binary systems will be in very close orbits, such that the red dwarf fills its Roche 
lobe. This leads to the transfer of mass, via an accretion disk, from the red dwarf to the 
white dwarf component - in other words, to cataclysmic variables (CVs). Yet searches for 
the many dozens of expected erupting CVs have typically found zero to two dwarf novae in 
clusters searched from the ground (Shara et al. 1994; 6 Tuairisg et al. 2003) and with the 
Hubble Space Telescope (HST) (Shara et al. 1996; Knigge et al. 2002). 

Recently, Chandra X-ray Observatory and HST observers claimed the detections of 
dozens of cataclysmic variable (CV) candidates in the cores of several globular clusters. 
These candidates are X-ray bright in Chandra observations (Grindlay et al. 2001) or 
UV-bright in HST observations (Knigge et al. 2002) and clustered in the cores of their 
host clusters. However, only zero or one of these UV or X-ray sources per cluster has yet 
exhibited the degree of variability (2 magnitude or larger outbursts) that originally led to 
labeling these variables as "cataclysmic." 

About half of all the CVs known in the solar neighborhood exhibit large amplitude 
outbursts — the Dwarf Novae. In the absence of spectroscopic confirmation (very 
challenging to obtain) and/or large amphtude outbursts, the X-ray and UV-bright CV 
candidates in globular clusters may remain just that: candidates, for the foreseeable 
future. It is thus important to fully substantiate every claim of large amplitude variability 
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in globular cluster stars. Only a few irrefutable, spectrographically confirmed CVs (e.g. 
Grindlay et al. (1995); Edmonds et al. (1999)), and cataclysmics confirmed by slitless 
far-UV spectroscopy (Knigge et al. 2003; Knigge 2004), (including 2 good candidates in 
Edmonds et al. (2003)) are known in globular clusters. A recent addition to this exclusive 
club of erupting CVs is the variable labelled VI in the dense globular cluster M15 (Charles 
et al. 2002, hereafter C02). 

Ml 5 is one of the densest globular clusters in the Galaxy. It has probably undergone 
core collapse and almost certainly contains a significant amount of non-luminous matter in 
its core (Baumgardt et al. 2003). The cluster hosts the remarkable low mass X-ray binary 
(LMXB) 4U2127-I-119, whose optical counterpart (Auriere et al. 1984) is the echpsing 
binary AC211 (Ilovaisky et al. 1987). M15 has been extensively imaged by HST, with 
which observers have discovered a sharp central density cusp (Guhathakurta et al. 1996), 
very blue stars and mass segregation in the core (De Marchi et al. 1996), Blue Stragglers 
(Yanny et al. 1994) and optically identified miUisecond pulsars (Anderson et al. 1990) — all 
indicators of strong dynamical interactions. If confirmed as an erupting CV, the candidate 
dwarf nova (C02) would be one of the rarest of exotic globular cluster stars, and an object 
worthy of further, intensive study. In anticipation fo our key result — that VI really is an 
erupting dwarf nova - we will hereafter refer to it as CVl. This is to avoid confusion with 
the nomenclature used in the Sawyer-Hogg (1973) catalog of variables in globular clusters, 
recently updated by Clement et al. (2001). There is an entirely different variable star 
labeled VI (in M15) in these catalogs. 

2. The M15 Dataset 

M15 has been imaged with the HST Wide Field Planetary Camera (WFPC2) during 6 
epochs, during one epoch with the HST Faint Object Camera (FOC), two epochs using the 
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HST Space Telescope Imaging Spectrograph (STIS), and finally with the HST Advanced 
Camera for Surveys (ACS) for two epochs. The wide range of filters and long time base 
of the observations (a decade) are particularly useful in characterizing variable stars. We 
downloaded all available images of M15 from the HST archive, and summarize the dataset 
in Table 1 and Figure 1. The variable is seen in eruption in the 4th and 5th epochs (see 
Figure 2), and in quiescence (below our frame limits) in the 1st, 2nd, 3rd, 6th, 7th, 8th, 9th, 
and the 11th epochs. It is detected in the near UV in the 10th epoch, possibly in transition 
between eruption and quiescence. There is Planetary Camera (0.044"/ pixel) coverage of 
CVl in the 1st, 2nd, 4th and ninth epochs. The 5th and 6th epochs have only WF2 (hence 
lower resolution, 0.1" /pixel) coverage. The 10th and 11th epochs were taken with ACS for 
our own survey for UV variables. Thumbnail images of the vicinity of CVl in most of the 
filters used in each epoch are presented in Figure 1; the variable is indicated in each frame. 



2.1. Photometry of the Dwarf Nova Candidate 

All of the WFPC2, ACS, and FOC data were obtained from the STScI Data Archive 
while the STIS images were still in the proprietary stage. All images were subject to the 
normal STScI on-the-fly recahbration pipehne software. The individual WFPC2 images 
were combined using the "Montage2" routine contained within the stand-alone DAOphot 
package. These frames were run through DAOphot's matching program "DAOmastcr" 
to derive the subpixel frame-to-frame shifts and thereby register all images to the first 
image of the first epoch. These shifts were passed to Montage2, which produced a single, 
sky subtracted, high signal-to-noise image for each epoch and filter. A variability search 
was performed by blinking rapidly through groups of four of the eleven epochs to look for 
any subtle changes in the appearance of the stellar field. Dwarf novae should rise from 
invisibility to easily detectable in the timeframe defined by the observations. Brightening 
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from previously empty regions of the sky is easily seen. Aside from CVl and the three other 
object mentioned in C02, no other DN-hke behavior was observed. 

Accurate measurement of the dwarf nova candidate's brightness is nontrivial because 
it is crowded between two neighboring stars. We have found that the most accurate way to 
carry out photometry on CVl is to use an aperture large enough to encompass CVl and its 
two companions. We measured the summed brightness of the two companions in epoch 1 
(when CVl is invisible), and verified in subsequent epochs (with CVl in quiescence) that 
these two stars do not vary. We subtracted off the brightness of these two companions in 
subsequent epochs. All of the results presented in this paper utilize aperture photometry 
while accounting for the varying platescales between instruments. 

The long-term light curve of CVl is shown in Figure 3. We have assumed a distance 
modulus to M15 of 15.0 (McNamara et al. 2004) and thereby transformed to the absolute 
magnitudes displayed in the Figure. The observations which use the ACS, FOC, and 
STIS instruments (with the F150LP, F165LP, F220W, F430W, and F28X50LP filters) are 
cahbrated to the ST magnitude system, while the WFPC2 magnitudes are in the Vega 
system. Vega zeropoints are only available for the WFPC2 filters, and so the ACS, FOC 
and STIS photometry was left in the ST system. Reddening corrections have been imposed 
on the photometry using a reddening value for M15 of E{B — V) — 0.09 (Cox 2000). Using 
this value, the reddening corrections in the Figs. 3 and 5 were A220 — 0.92, Au — 0.41, 
and Ab = 0.36 for the F220W, F336W, and F439W filters, respectively. Similar reddening 
corrections, appropriate for each band, were applied to the derived platelimits in Fig. 
3. CVl was fainter than the WFPC2 frame limit corresponding to absolute magnitude 
~ 6.0 in each of the April 1994 and August 1994 images. The variable reached maximum 
brightness as seen in the F336W filter on 26 October 1994, at absolute magnitude +1.5. It 
was next imaged in December 1998 for seven consecutive days in multiple filters, when it 



-8- 



was also in eruption. We have assembled the photometry for every available image during 
this week in Figure 4. There are indications that the variable was dimming towards the end 
of this epoch. CVl was below our detection threshold in the remaining epochs, except for a 
detection in the near UV with the F220W filter on ACS on 27 October 2003. 

3. The Eruptions of CVl 

Charles et al (2002) discovered CVl in eruption in U band (F336W) HST PC images 
of October 1994 (seen in Figure 1 as epoch 4). They predicted that the object is blue, and 
suggested (since X-ray transients erupt only rarely) that CVl is much more likely to be a 
dwarf nova than an X-ray transient. 

The December 1998 images (epoch 5) all have CVl on the (lower resolution) HST 
WF. CVl is somewhat blended with its nearest neighbors, but careful measurement of the 
position of CVl (see also Figure 2) shows that it is resolved and seen in eruption in F336W 
during the 1998 observations (contrary to the claim by C02). Comparison of the F336W 
and F439W images of CVl in this epoch (Figure 1) clearly show that CVl is very blue, 
with U — B —0.53 as predicted by C02. The hght curve (and peak eruption luminosity) 
of CVl shown in Figures 3 and 4 are very similar to the typical eruption light curve and 
luminosity of the prototypical dwarf nova SS Cygni centered on Julian Date 24501226, 
shown in Figure 5. The U — B color is also similar to that of SS Cyg (Hopp & Wolk 1984) 
in outburst. The apparently significant fiickering seen in the F502N and the Ha images 
(Figure 4) is often seen in erupting dwarf novae. 

SS Cygni erupts, on average, every 40 days, and is brighter than minimum light for 16 
days during this time (Szkody & Mattel 1984). It is nearly at maximum for 7 days, followed 
by a 7 day long decline to minimum. If CVl is like SS Cygni then we should see it above 
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its minimum state ~ 16/40 = 40% of the time. In fact we detect CVl in eruption or above 
minimum in 3 of 11 epochs, or 27% of the time — an entirely consistent result, given the 
small number statistics. 

Finally, CVl is faint but visible in the F220W images of epoch 10, at absolute 
magnitude 5.7, about 4.0 magnitudes fainter than in epochs 4 and 5 (ignoring color effects). 
If we assume a UV color (F220W - C/) ~ —0.5 then CVl is still 3.5 magnitudes fainter in 
epoch 10 than in epochs 4 and 5. We suggest that this is an intermediate state between the 
star's quiescent and outburst states. This is because CVl is NOT seen in the far-UV images 
(of epoch 11), which should be even more sensitive to UV-bright CVs than the near-UV 
imagery of epoch 10. 



3.1. Is CVl a Dwarf Nova? 

The absolute magnitude of CVl at maximum light, coupled with its very blue color, 
similarity in light curve to SS Cygni near maximum, detection in near UV but not in visible 
light near minimum, relative frequency of eruption (at least twice and possibly three times 
in eleven independent epochs) and amplitude of eruption (at least 4 and possibly more 
magnitudes) are all strongly supportive of CVl being a dwarf nova. 

While the case is not perfect (we do not have confirmatory spectra or a full light curve 
from quiescence to outburst and back again) it is very strong, and we fully concur with the 
characterization of C02 of CVl as a dwarf nova. 

We predict that an HST outburst spectrum of CVl would show Balmer absorption 
lines, but both the crowding of CVl in particular and the unpredictability of dwarf nova 
outbursts in general will make this a challenging observation (the S/N for the the Ha 
imagery of epoch 5 is too low to test this prediction). 
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While Ml 5 now appears to possess a dwarf nova, it remains the only one known in this 
dense, massive cluster. The two week duration search for dwarf novae in M15 by O Tuairisg 
et al. (2003) failed to detect any new candidates. The remarkable result that either zero, 
one or at most two erupting dwarf nova can be found in massive globular clusters seems to 
be holding true. 
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Fig. 1. — A matrix of HST images of the field of the candidate dwarf nova CVl in M15. 
North and East are as indicated, and each postage stamp image is 1.24 arcsec on a side. The 
upper labels are the HST filter names. CVl is indicated in each postage stamp image. 
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Fig. 2. — A contour plot showing the neighborhood of the erupting Dwarf Nova CVl. From 
left to right, the plots show the first (quiescent), and fourth and fifth (eruptive) epochs, 
respectively in F336W. The arrow shows the location of the Dwarf Nova. Each plot is 1.09 
arcsec on a side. The contours of the left and middle plots indicate 40-count intervals, while 
the right plot shows 20 counts per contour. 
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Fig. 3. — Eleven epochs of multi-wavelength photometry of CVl in M15. The object is in 
eruption in epochs 4 and 5, and probably in an intermediate state in epoch 10. Also shown 
are the upper-detection limits for those epochs in which CVl is quiescent. 
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Fig. 4. — Daily observations of the dwarf nova CVl in outburst over 1 week in December 
1998 (epoch 5). Magnitude offsets are arbitrary. 
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Fig. 5. — The dwarf nova SS Cygni in outburst over two different time resolutions (upper 
two plots). Observations are visual, from the AAVSO. Also shown is the lightcurvc of CVl 
in M15 (same as Figure 4) in F336W and F439W. The light curves suggest that CVl was 
in early decline from a dwarf nova outburst in epoch 5. 



Table 1. Table of M15 observations: Listed are the number of useable observations at each date. The column labels 

from F150LP to F814W are HST filter names, in order of increasing central wavelength. Also, for the WFPC2 
images, the placement of the center of the cluster falls either on the PC chip or a WF, as indicated. 
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